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Abstract 

A Fortran program for the calculation of the ground state properties of axially 
deformed even-even nuclei in the relativistic framework is presented. In this rela- 
tivistic mean field (RMF) approach a set of coupled differential equations namely 
the Dirac equation with potential terms for the nucleons and the Glein-Gordon type 
equations with sources for the meson and the electromagnetic fields are to be solved 
self-consistently. The well tested basis expansion method is used for this purpose. 
Accordingly a set of harmonic oscillator basis generated by an axially deformed 
potential are used in the expansion. The solution gives the nucleon spinors, the 
fields and level occupancies, which are used in the calculation of the ground state 
properties. 



1 PROGRAM SUMMARY 



Title of program: RMFD.f 
Catalogue number: ?? 
Program available from: 
Licensing provisions: none 
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Computer for which the program is desinged and others on which it has been tested: any 
Unix work-station or mainframe or PC. The program has been tested on work stations: 
DEC, DEC- Alpha , on CDC6600 mainframe and on 486 IBM compatible PC's 

Operating system: UNIX or VMS or MS-DOS 

Programming language used: Fortran 77 

No. of lines in distributed program, including test data, etc: ca. ?? 

Keywords: Relativistic mean field theory, binding energy, nuclear radii, deformations and 
densities. 



Nature of physical problem 

The Relativistic Mean Field (RMF) theory [I], 0] has been astonishingly succesfull || [|, 
|5], |6], [7|, H in accurately describing the nuclear matter properties and the ground state 
properties of finite nuclei spread over the entire periodic table including those away from 
the stability line. The nucleonic and mesonic degrees of freedom are explicitly included 
from the very beginning in the relativistic framework. As a result the correct spin-orbit 
splittings emerge automatically. Initially the RMF theory was applied succesfully to the 
description of the properties of spherical nuclei. As most of the nuclei are deformed, the 
generalization of the solution of the RMF equations for this case is required, which is 
a non trivial task. Therefore a computer program was developed ||, 10, H to solve the 
RMF equations, suitable for the calculation of the ground state properties of the axially 
deformed nuclei. The present program is an improved version and also is compatible for 
PC's. 



Method of solution 

In the RMF theory one needs to solve self-consistently a set of coupled equations namely 
the Dirac equation with potential terms for the nucleons and the Klein-Gordon type equa- 
tions with sources for the mesons and the photon. For this purpose we employ the well 
tested basis expansion method. The bases used here, are generated by an anisotropic 
(axially symmetric) harmonic oscillator potential. The upper and lower components of 
the nucleon spinors, the fields as well as the baryon currents and densities are expanded 
separately in these bases. The expansion is truncated so as to include all the config- 
urations up to a certain finite value of the major oscillator shell quantum number. In 
this expansion method the solution of the Dirac equation gets reduced to a symmetric 
matrix diagonalization problem, while that of the Klein Gordon equation reduces to a 
set of inhomogenous equations. The solution provides the spinor, fields, and the nucleon 
currents and densities (sources of the fields), from which all the relevant ground state 
nuclear properties are calculated. 

Restrictions on the complexity of the problem 

The present version is applicable to even-even nuclei due to the imposition of time reversal 
inavariance and charge conservation. The program can be modified for the general case 
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including that of odd mass nuclei by incorporating the additional currents arising due to 
time reversal breaking. 

Typical running time: 

From 30 minutes to several hours depending upon the computer for the general case. 
However, the computer time will considerably increase if one wishes to include higher 
number of shells. 

2 LONG WRITE-UP 

2.1 RMF equations 

The basic Ansatz of the RMF theory is a Lagrangian density |T|, Q where nucleons are 
described as Dirac particles which interact via the exchange of various mesons. The 
Lagrangian density considered is written in the form: 

C = ^(ip - M)i/j + ±<V^(x - U(a) - \n ia ,sv w + 

\mlu^ - IR^BT + \m 2 p p^ - \F^F^ (1) 

ga^crip - g^ipip - g P ^pr^ - eipflip 

The meson fields included are the isoscalar a meson, the isoscalar-vector u meson and 
the isovector-vector p meson. The latter provides the necessary isospin asymmetry. 

The arrows in Eq. (1) denote the isovector quantities. The Lagrangian contains also 
a non-linear scalar self-interaction of the a meson. 

U(a) = -mla 2 + ~g 2 a 3 + \g,a A (2) 

This term is important for appropriate description of surface properties IjTTJ. M, m CT , m^ 
and nip are the nucleon-, the a-, the uj- and the p-meson masses respectively, while g CT , 
g w , g p and e 2 /47r = 1/137 are the corresponding coupling constants for the mesons and 
the photon. The field tensors of the vector mesons and of the electromagnetic fields take 
the following form: 

W = d»uo u - d v u» (3) 
W v = - d v p» (4) 

The variational principle gives the equations of motion. The mean field approximation 
is introduced at this stage by treating the fields as the c-number or classical fields. This 
results into a set of coupled equations namely the Dirac equation with potential terms for 
the nucleons and the Klein-Gordon type equations with sources for the mesons and the 
photon. For the static case, along with the time reversal invariance and charge conserva- 
tion the equations get simplified. The resulting equations, known as RMF equations have 
the following form. 
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The Dirac equation for the nucleoli: 



{-iaV + V(v) + fi[M + S(v)}}^ = e^i, (6) 
where V(r) represents the vector potential: 

1 + r 3 

V(r) = ^cuo(r) + g p r 3 p (r) + e^— A (r), (7) 
and S(r) is the scalar potential: 

S(r) = gM?) (8) 
the latter contributes to the effective mass as: 

M*(r) = M + S(r). (9) 

The Klein-Gordon equations for the meson and the electromagnetic fields with the 
nucleon densities as sources: 

{-A + m>(r) = -g aPs (r) - g 2 a 2 (r) - g 3 a 3 (r) (10) 

{-A + m^ (r) = g uPv (r) (11) 

{-A + mJ}p (r) = g p p 3 (r) (12) 

-AA (r) = ep c (r) (13) 

The corresponding densities are: 



Ps = 


A 

■i=i 




Pv = 


A 

i=i 




P3 = 


z 

E nii>+ ip p - 
p=l 


N 

E i) n 

n=l 


Pc = 


z 

E niip+ ifj p . 
P =i 





(14) 



Here the sums are taken over the particle states only. This implies that the contributions 
from negative-energy states are neglected (no-sea approximation), i.e. the vacuum is not 
polarized. The n meson does not contribute in the present relativistic mean field (Hartree) 
approxiamtion because of its pseudo nature. The occupation number n« is introduced to 
account for pairing which is important for open shell nuclei. In the absence of pairing it 
takes the value one (zero) for the levels below (above) the Fermi surface. In the presence 
of pairing the partial occupancies (rij) are obtained in the constant gap approximation 
(BCS) through the well known expression: 

n t = 1(1 ; El ~ A ) (15) 



4 



The £j is the single-paricle energy for the state i and chemical potential or Fermi energy 
A for protons (neutrons) is obtained from the requirement 



rii = the number of protons (Z) 

(16) 

(the number of neutrons (N)) 

The sum is taken over protons (neutrons) states. The gap parameter A is calculated from 
the observed odd-even mass differences. In the absence of experimental masses it can 
inferred from the extrapolation of the masses given by any of the avalaible mass formulae 
(e.g from ref. ||12|| ) 



The above set of equations (6,10-13) are to be solved self-consistently. For this purpose 
one starts with an initial guess of the fields (e.g. generated by axially deformed Woods- 
Saxon potential) to calculate the potential terms (Eqs. (7,8)) appearing in the Dirac 
equation (Eq. (6)). The Dirac equation is solved with these potentials terms to yield 
the nucleon spinors which in term are used to obtain the sources (densities). The meson 
and photon equations are then solved with these sources to get a new set of fields to be 
used for the calculation of new potential terms. The Dirac equation is then solved with 
the new potentials to get the spinors again to be used to obtain the new sources for the 
meson fields. This iterative procedure is continued till the converegence upto the desired 
accuracy is achieved. 



2.2 Axially symmetric case 

For the axially symmetric deformed shape the rotational symmetry is broken and there- 
fore, the total angular momentum j is no longer a good quantum number. However, the 
densities are still invariant with respect to a rotation around the symmetry axis, which 
is taken to be the z-axis. It then turns out to be useful to work with the cylindrical 
coordinates 



x = r ± cos (p, y = r ± sin (p 



and 



(17) 



For such nuclei the Dirac equation can be reduced to a coupled set of partial differential 
equations in the two variables z and r x . In particular, the spinor ipi with the index i is 
now characterized by the quantum numbers 



5 ^ii ti 



where is the eigenvalue of the symmetry operator j Zi (the projection of ji on the z-axis), 
7Ti is the parity and ti is the isospin. The spinor can be written in the form: 



(ft 

fi 

ig'i 



-l/2)<p v 



xtM 



(19) 
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The four components fi(r ± ,z) and g^(r ± ,z) obey the Dirac equations 



(M* + V)f+ 


+ d z gt 


+ (d r± 


Q + 
+ 

r ± 


i 


= 


(20) 


(M* + V)fr 


- 9 z Qi 


+ (d r± 


Q- 

r ± 


-W 


= e% fr 


(21) 


(M* - V)gt 


+ d z ft 


+ (dr ± 


n + 
+ — 

r ± 






(22) 


(M* - V)gr 


- djf 


+ (9r ± 


Q - 
r, 




= -£i 


(23) 



For each solution with positive Q: 

V* ={f?,fr,g?,gr,ni} (24) 

we have the time reversed solution with the same energy 

^ = = {-frj^g-^gt, -OJ (25) 

with the time reversal operator T = zcr y i^ (if being the complex conjugation). For nuclei 
with time reversal symmetry, the contributions to the densities of the two time reversed 
states i and i are identical. Therefore, we find the densities 

Ps , v = 2j>((i ft r + 1 fr i 2 ) t (i & + r + 1 ft- 1 2 )) (26) 

and, in a similar way, p% and p c . The sum i > runs only over the states with positive 
Qi— values. These densities serve as sources for the fields = a, u° p° and A , which are 
determined by the Klein-Gordon equation in cylindrical coordinates: 

( — -d r± r ± d r± - d 2 + m/)4>(z,r ± ) = s^{z,r ± ) (27) 
fx 

The inhomogeneous parts are given by 

-g a p s {z, r ± ) - g 2 cr 2 (z, r ± ) - g 3 cr 3 (z, r ± ) for the tr-field 

r , = ^ gcPv(z,r ± ) for the ^-field ^ 

" g p p 3 (z,r ± ) for the p— field 

ep p (z,r ± ) for the Coulomb field 

2.3 Solution of the RMF equations 

For the solution of thr RMF equations we use the basis expansion method. We follow 
closely the details, presentation and the notation of ref. ||. For the axially symmetric 
case we expand the spinors f i and gf in Eqs. (20-23) in terms of the eigenf unctions of a 
deformed axially symmetric oscillator potential 

V osc {z,r ± ) = -Mu 2 z z 2 + -Mooy ± (29) 



Imposing volume conservation, the two oscillator frequencies hu ± and huo z can be ex- 
pressed in terms of a deformation parameter (3$: 

hu z = hu Q exp(-y^/3 ) (30) 
hu ± = huj exp(+±y/£(3 ) (31) 

The corresponding oscillator length parameters are 



h = K = (32) 

Because of volume conservation, we have b 2 ± b z = b^. The basis is now determined by 
the two constants Tiojq and The eigenfunctions of the deformed harmonic oscillator 
potential are characterized by the quantum numbers 

| a > = | n z ,n r ,mi,m s > (33) 

where mi and m s are the components of the orbital angular momentum and of the spin 
along the symmetry axis. The eigenvalue of j z , which is a conserved quantity in these 
calculations, is 

Q = mi + m s- (34) 

The parity is given by 

vr = (35) 
The eigenfunctions of the deformed harmonic oscillator can be written explicitly as 

$ a (z,r ± ,<p,s,t) = <p nz (z) ^:( r ^ eimiiP Xms(s)Xt a (t) (36) 



with 



*a(r,s)xt a (t) (37) 



fjrni 

CKrx) = -f- V2n m ^L2( V )e-^ (39) 



and 



C = z/b x , r, = r 2 Jbl (40) 

The polynomials H n (() and L™{rf) are Hermite polynomials and associated Laguerre 
polynomials as defined in ref. [13[ The normalization constants are given by 



1 n I 

N nz = , = and AC = J , f r. (41) 



In order to evaluate the matrix elements, we also need the polynomials L™ 1 and H Uz 
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defined by the derivatives 

d z <PnXz) = §f^(C)e~ c2/2 (42) 
bz 

N m i 

drJQ(r x ) = -^V^n^-^L^e-^ (43) 

which can be calculated from the recursion relations 

H n XQ = (H nz (()-H nz+1 (() (44) 
L%(n) = (2n r + m l -r ] ) L%(n) - 2(n r + m,) L^-M (45) 

The solutions of the Dirac equation in the axially symmetric case have only the good 
quantum numbers Q and 7r and we use the expansion 

fi(*,*,t) = ^(^I;j^J = E f^^a(r,s) X tM (46) 
9i{*,s,t) = -T^l 9 !- Y:l !i(n+i/2)J = E 9$ ^(r,s) X tM (47) 



To avoid the appearence of the spurious states || [14| the quantum numbers a max and 
aWr are chosen in such a way that the corresponding major quantum numbers N = 
n z + 2n p + mi are not larger than Np + 1 for the expansion of the small components, 
and not larger than Np for the expansion of the large components. The Dirac equation 

reduces to a symmetric matrix diagonalization problem: 

■£;' %)(%) - «.(!) <«> 

of the dimension ct max . + a max . The matrix elements ^4 aa ', B aa > and C aQ / are given by: 

•Acta 1 
C aa i 



5 mm >5m s m' s NZ l N nz N™\N n , J d V e-^ m 'L^( V )LZ(v) (49) 

o 

Jd(e-( 2 H nz (()H n A0 (M*(b z (,b ±v ^j) ± K(6zC,6-lV^)) (50) 



oo 

X 





n -x a a (-)^" ms ^ i^ 7 A I , 

Oqo' — mimi i0 msrri i g nrnr i \°n z 'n z +l]l ~^ °n z n z '+iy ] l'" ) 

N n r ' N Z 1 ' f 

(52) 

oo 

/ d V e-^ m ^ 2 L^(n) ~ m^iv)) (53) 

o 

<W.+i ]dr,e-^ mi - 1/2 L7M + (mi + !)!#(»?))} (54) 



V s m s +1 
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In the next step, we calculate the density matrices in the shell model space: 



p (i) f (i) n (i) n (t) 

i>0 

and in coordinate space: 



p'£ = zEmifM =f (55) 

i>0 



ps,v( z , r ±) = ^ " Z & (56) 

x N nz N n ,H n XQH n ,{<) K;N™\rT l L^{ V )L^{v)- (57) 

Similarly, one obtains p3 and p c . These are the sources for the solution of the Klein-Gordon 
equations. The fields of massive mesons are also solved by the expansion in a deformed 
oscillator basis. Here for computational and numerical convenience we use the same 
deformation parameter j3 as in Eqs (32,33), and take the oscillator length = bo/y/2 : 



N - 

(:./•_ 

X >B n z n r 



1 " S 

e -C 2 /2-,/2 £ <P nzn N nz H nz (C) V2Ll( V ) (5f 

Of 



with ( = \/2z/b z and rj = 2r\jb\. Inserting this ansatz into the Klein-Gordon equation, 
we obtain an inhomogeneous set of linear equations 

N B 

?j 7~Ln z n r n z 'n r ' ( t ) n z 'n r ' = S n z n r (59) 
n z 'n r ' 

with the matrix elements 

2 2 

H nznrnz i nr i = (- — (n z + l/2) - —(2n r + l) + m 2 ^)5 nrnr >5 nznzl (60) 
+ rkV in x + l)nj 5 nz > nz+1 + -J(n z ' + l)n z 5 nzriz > +1 ) (61) 

b z 2 v 2 v 

2 

_ Tr( n r <Wn r +l + n r 5 nrnr > +1 ) (62) 

For the Coulomb field, due to its long range character we cannot use the oscillator basis 
expansion method. Therefore for this case the standard Green's function method is used 



as is discussed in the appendix of ref. |L5 



In the present case with time reversal symmetry and pairing the total binding energy 
is given by the sum of various individual contributions: 

E= E part + E a + E u , + E p + E c + E pair + E CM -AM (63) 

with 

E P art = ^2ni 2 ei (64) 

i 

E aL = "f/ ^p s (r) ff (r) (65) 



Knl = ~\ J 'd 3 r {^ 2 a(if + ^(r) 4 } (66) 

E u = "f / d 3 rp v (v)cu°(r) (67) 

E p = -|/rf 3 rp 3 (r)p 00 (r) (68) 

2 

E c = -|- / rf 3 rp c (r)A°(r) (69) 

£ poir = -A ^ \/ni(l - 77-i) (70) 

i 

3 3 i 

E C m = --hio = --4L4-3 (71) 



In the calculation of the pairing energy we use a pairing window, i.e. the sum over i in 
Eq. (69) is only extended up to the level where — A < 2(41v4~ »). The factor 2 has been 
determined so as to reproduce the pairing correlation energy (Eq. (69)) for neutrons in 
the nucleus 118 Sn calculated by using the Gogny force [PSfl . 

The charge radius is calculated using the following formula: 

r c = + 0.64 (fm) (72) 

The factor 0.64 in Eq. (71) accounts for the finite size effects of the proton. 

The quadrupole Q n , v and hexadecupole H np moments for neutrons and protons are 
calculated using the expressions: 

Q n , p = < 2r 2 P 2 (cos6) > n , p = < 2z 2 - x 2 - y 2 > n , p (73) 

and 



H n>p = < r 4 Y 40 (9) > n , p = y — < 8z 4 - 2Az 2 (x 2 + y 2 ) + 3(x 2 + y 2 ) 2 > n , p (74) 

The conventional deformation parameter (3 is obtained from the calculated quadrupole 
moments through 

F\(hr 3 

Q = Q n + Q P = ^——AR 2 Q (3. (75) 

with Rq = 1.2A5 (fm). 

2.4 Calculation of ground state properties 

The present program has been used to calculate the ground propetries of axially deformed 
even-even nuclei spread over the entire periodic table encompasing also the nuclei far away 
from the stability line. 

For carrying out explicit numerical calculations one requires the following input infor- 
mation: 
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i) The baryon and the meson masses and the coupling constants of the meson fields 
to the nucleons appearing in the Lagrangian (Eq.(l)) 

ii) The number of oscillator shells Np and Nb, i.e. the cutoff major oscillator shell 
quantum number up to which the Dirac spinors (Fermion wavefunctions) and meson fields 
(describing the Bosonic degrees of freedom), and so also the densities, are expanded. 

iii) The basis parameters hu and /3q used for the expansion of the Dirac spinors and 
the fields. Following ref. || we fix for fermions hto = 41v4~ 1//3 and (3q can be taken to be 
any reasonable value (preferably close to the experimental value). 

As stated before the oscillator size parameter for the meson fields (and the densities) 
are fixed at l/v2 times the corresponding Fermion parameters. The deformation pa- 
rameter of the oscillator basis (3q has been chosen to be identical for the Dirac spinors 
and the meson fields. This simplifies the calculations and avoids the need for additional 
parameters. 

For illustration we present here the results of calculations for the ground state prop- 
erties of Sr nuclei over a wide range of isospin [0, namely Sr isotopes with mass numbers 
A=70 up to A=110. The calculations have been carried out using the Lagrangian param- 



eter set NL-SH, which provides good results on both sides of the stability line |I7| . The 
NL-SH set of parameters are listed in Table 1. 

The number of shells taken into account in the expansion is 12 for both fermions and 
bosons (Nf=Nb=12). It should be noted that for convergence reasons 14 shells were also 
considered. It was observed that, there is hardly any difference between these two sets of 
results. This observation also holds for nuclei in the entired periodic table. Therefore, it 
is sufficient to consider 12 shells in practical calculations. However, for superheavy nuclei 
probably one has to include higher number of shells. 

For open shell nuclei, pairing has been included using the BCS formalism. In the 
BCS calculations we have used constant pairing gaps, which are taken from the empirical 
particle separation energies [18|] of neighbouring nuclei. 



Now we present and discuss some of the calculated results. Fig. 1 shows the binding en- 
ergy per nucleon (E/A) for Sr isotopes. The empirical values taken from the 1993 Atomic 
Mass Evaluation Tables fl8| (expt.) are also shown. The figure also includes the predic- 
tions of the recent finite-range droplet model |L2| (FRDM) and of the extended Thomas- 
Fermi with Strutinsky Integral (ETF-SI) model |HJ for comparison. The parabolic shapes 
of the binding energy per nucleon emerges nicely. The minimum in the binding energy 
is observed at the magic neutron number N=50 in the RMF as well as in the mass mod- 
els. The calculated RMF binding energies agree closely to the predictions of the FRDM 
and ETFSI models within 1-2 MeV. The RMF theory predicts binding energies which 
are in accord with the empirical values in almost all the cases, with deviations < 0.5%. 
It should be noted that similar type of agreement has also been obtained in other mass 
regions. 

We have performed calculations in the RMF theory for both the prolate and oblate con- 
figurations. The deformations of nuclei have been obtained from the relativistic Hartree 
minimization. We show in Fig. 2 the quadrupole deformation (3 2 for the shape corre- 
sponding to the lowest energy. The predictions of FRDM and ETF-SI are also shown for 
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comparison. It is seen that the RMF theory gives a well-defined prolate shape for lighter 
isotopes. Further, an addition of a few neutrons below the closed neutron shell leads to an 
oblate shape. This shape turns into spherical ones as nuclei approach the magic neutron 
number N=50. Nuclei above this magic number revert again to the prolate shape in the 
RMF theory. Thus, a shape transition from prolate-oblate-spherical-prolate is followed. 

In addition to the lowest minimum, several isotopes exhibit a second minimum, thus 
implying a shape-coexistence, i.e., the prolate and the oblate shapes differ in the energy 
only by a few hundred keV. These nuclei have been shown by squares surrounding the 
black circles. The calculations predict two minima for several heavy Sr nuclei, the prolate 
shape results being a few hundreds keV lower in energy than the oblate one. This is 
displayed in Fig. 3, where the difference in the ground-state binding energy of the oblate 
and the prolate configurations is shown. It can be noticed that Sr isotopes beginning 
with A=92 acquire a prolate shape predominantly. For nuclei close to A=98, the prolate 
shape is lower than the oblate one only by about 300 keV. With a further increase in the 
neutron number the Sr isotopes take up the prolate minimum, the oblate shape being 
about 600-800 keV higher. The undulation in the prolate-oblate energy differences of 
neutron-rich Sr isotopes is a noteworthy feature of the RMF prediction. 

The (3 2 values from NL-SH, FRDM and ETF-SI is shown in Table 2. The (3 2 values 
for the second minimum obtained for Sr isotopes with the force NL-SH are shown in the 
parentheses. A comparison shows that the three approaches provide values close to each 
other. The experimental quadrupole deformations obtained from BE(2) values taken from 



22fl are shown in the last columns of the table. It may be noted that these empirical /?2 
values do not indicate any signature as to ascertain the shape of a given nucleus. The 
absolute values, however, do compare well with the RMF predictions. 

In Fig. 4 the rms charge and neutron radii of Sr nuclei are presented. It is seen that on 
going from the lighter isotopes to the heavier ones the charge radii exhibit a decreasing 
trend upto the magic isotope, that is the lighter isotopes have higher charge radii than the 
heavier closed neutron- shell nucleus. The charge radii for nuclei heavier than the closed 
neutron-shell start increasing with addition of neutrons. The neutron radii, on the other 
hand, also show a kink about the neutron shell closure. However, the neutron radius for 
lighter isotopes in these chains is not higher than that of the closed-shell nucleus. 

In Fig. 5 we show the isotope shifts r^(A) — r^(re/)) for Sr nuclei calculated with 
respect to a reference nucleus ( 88 Sr). The empirical data obtained from atomic laser 
spectroscopy |20, |21| are also shown. The experimental data for Sr nuclei exhibit a kink 
about the magic neutron number. This kink about the closed-shell is a characteristic 
feature of isotope shifts in many nuclei P0|. A solution to this problem has eluded since 



long. It can be seen from the figure that the RMF theory is successful in reproducing 
this kink. We would like to stress that the RMF theory provides a first-ever microscopic 
description of anomalous isotopic shifts in Sr nuclei. Such an anomalous behaviour is a 
generic feature of deformed nuclei which include almost all isotopic chains in the rare- 
earth region [f20| . It is to be remarked that here also the RMF theory is shown to have a 
remarkable success ||. This is amply demonstrated through Fig.6 where the isotope-shift 
of rare-earth nuclei with atomic numbers ranging from Z=60 upto Z=70 are shown. 
In short the RMF theory is highly successful in describing the ground state properties 
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( e.g binding energies, nuclear radii, the deformations, isotope shifts etc). The agreement 
with the experiment is really remarkable. It is to be stressed that the RMF theory uses 
limited number of parameters as compared to other theories and models and at the same 
times it provides a consistent and unified description of the ground state properties of 
nuclei over range of isospin. 

3 Program structure and test run 

The code consists of the fortran program and two additional files: DIZ.PAR and DIZ.DAT. 
The file DIZ.PAR contains the relevant information regarding the dimensions, depending 
upon the number of oscillator shells to be included. The standard form of DIZ.PAR is for 
12 shells both for fermions and bosons (Nf—Nb—12). The input file DIZ.DAT provides 
the necessary information relevant to the specific case being calculated (the details are 
attached). The program runs interactively. After three iterations the user has to provide 
as input the number of iterations and the value of the parameter xmix. As an example for 
hundred iterations and the value of xmix=0.3 the user should write -100 and in the next 
0.3. The main program DIZ calls various subroutines reading the data and performing 
the execution. The operation consists essentially of two parts. The first part using 
DIZ.DAT starts the program, initializes and generates all the relevant information. It uses 
the subroutines PREP, READER, START, DEFAULT, GAUSS. The fortran file DIZLIB 
contains several general subroutines required at various stages in the program. The second 
part is the main part and does the entire execution using the initial information provided 
by the first part. The iterative procedure is carried out by the subroutine ITER. In the 
first iteration it solves the Dirac equation using the potential terms calculated by the 
initial guessed values of the fields. The solutions (Dirac spinors) are used to calculate the 
sources which in turn are used in solving the Klein Gordon equations. These solutions 
are used in the next iteration for the solution of the Dirac equation. This procedure 
continues till the convergence of the desired accuracy is obtained. In this operation ITER 
calls various subroutines like POTGH, DIRAC, FIELD, OCCUP, DENSIT, EXPECT 
etc. The output file (DIZ. OUT) is prepared in RESU and INOUT. The user is advised 
to rerun the program once more after getting the convergence. This will make the output 
file DIZ. OUT compact. The file DIZ. OUT contains explicit headings to make it self- 
explainatory. Several detailed comments are introduced at various places in the program 
which helps the user to understand the different functions of subroutines and at the same 
time to figure out what is going on at important steps of the program. 

For an illustration DIZ.DAT file is listed for the specific case of 88 Sr. A part of the 
output file DIZ. OUT is also listed providing the results of the RMF calculation for this 
nucleus, using the parameter set NL-SH. If ones wishes to use a different set of Lagrangian 
parameters then the new set of parameters can be inserted in place of the corresponding 
numbers of the NL-SH set in the subroutine DEFAULT. 
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Figure Captions. 



Fig. 1 The calculated RMF binding energy per nucleon for Sr isotopes obtained with 
the parameter set NL-SH. The predictions from the mass models FRDM and ETF-SI are 
also shown for comparison. 

Fig. 2 The quadrupole deformation fa obtained from relativistic Hartree minimization 
for Sr isotopes using the force NL-SH. The predictions of the mass models FRDM and 
ETF-SI are also displayed for comparison. Nuclei exhibiting a shape coexistence and thus 
a second minimum in the RMF theory are depicted by a square surrounding the fa value 
of the lowest minimum. 

Fig. 3 The prolate-oblate shape coexistence for neutron-rich Sr isotopes predicted in 
the RMF theory. The energy difference in the prolate and oblate minima for Sr isotopes 
is shown. 

Fig. 4 The calculated rms charge and neutron radii of Sr isotopic chain obtained by 
using the force NL-SH. 



Fig. 5 The calculated and the experimental |20|, |21] isotope shifts for Sr isotopes. 



Fig. 6 The calculated and the experimental isotope shifts for Nd, Sm, Gd, Dy, Er, Yb 



nuclei. The empirical values for all but Gd nuclei are taken from ref. pOfl. The isotope 



shifts for Gd nuclei have been derived from ref. 123 . 
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Table 1: The parameters of the force NL-SH. All the masses are in MeV, while gi is in fm 1 . 
The other coupling constants are dimensionless. 



M = 


939.0 


m c = 526.059 




= 783.0 


m p = 763.0 


9a = 


10.444 


9u, = 12.945 


9 P 


= 4.383 


g 2 = -6.9099 £ 3 = -15.8337 
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Table 2: The quadrupole deformations /3 2 for Sr isotopes obtained in the RMF theory using 
the force NL-SH. The FRDM and ETF-SI predictions are also listed. The available empirical 
deformations (expt.) obtained from the BE(2) values are also given in the last column. The 
experimental values do not depict the signature of the deformation. The deformations of 
nuclei showing a shape coexistence with a second minimum are given in the parentheses. 



A 


NL-SH 


FRDM 


ETF-SI 


expt. 


72 


0.324 


0.371 


-0.30 




74 


0.430 


0.400 


0.44 




76 


0.450 




0.421 


0.44 - 


78 


0.450 


0.421 


0.43 


0.434 


80 


0.402 


0.053 


0.40 


0.377 


82 


-0.200 


0.053 


-0.30 


0.290 


84 


0.089 


0.053 


0.15 


0.211 


86 


0.0 


0.053 


0.00 


0.128 


88 


0.0 


0.045 


0.00 


0.117 


90 


-0.058 


0.045 


-0.11 




92 


0.181(-0.165) 


0.080 


-0.15 




94 


0.230(-0.218) 


0.255 


-0.19 




96 


0.356(-0.275) 


0.338 


0.35 




98 


0.424(-0.309) 


0.357 


0.39 


0.354 


100 


0.426(-0.314) 


0.368 


0.38 


0.372 


102 


0.413(-0.295) 


0.369 


0.40 




104 


0.403(-0.277) 


0.361 


0.38 
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